Collapse of Bose-Einstein condensate with dipole-dipole interactions 
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A dynamics of Bose-Einstein condensate of a gas of bosonic particles with long-range dipole- 
dipole interactions in a harmonic trap is studied. Sufficient analytical criteria are found both for 
catastrophic collapse of Bose-Einstein condensate and for long-time condensate existence. Analytical 
criteria are compared with variational analysis. 



PACS numbers: 03.75.Fi, 05.30.Jp 



(N 
O 
O 
(N 

< 

in 



> 

m 
oo 
O 
n 
o 



o 
o 



Bose-Einstein condensation of dilute trapped atomic 
gases [|[ p| essentially depends on the interpaticle inter- 
actions. In most experiments so far the dominated inter- 
actions were short-range van der Waals forces which are 
characterized by the s— wave scattering length a. Spa- 
tially homogeneous condensates with positive scattering 
length (repuilsive interaction) are stable while conden- 
sates with negative scattering length (attractive interac- 
tion) are always unstable to local collapses (|] because 
the quantum pressure is absent in homogeneous conden- 
sates. The presence of trapping field allows to achieve a 
metastable Bose-Einstein condensate (BEC) Q for a < 
if the number of particles is small enough to ensure exis- 
tence of local minima of energy functional [[| . 

Recent progress in creating of ultra-cold molecular 
clouds |, §] opens a new prospective for achieving BEC 
in a dilute gas of polar molecules and stimulates growing 
interest in study of BEC with dipole-dipole interactions 
d, |, |, HI |l|. Dipole-dipole forces are long- 
range and essentially anisotropic. Net contribution of 
dipole-dipole interactions can be either repulsive (posi- 
tive dipole-dipole interaction energy) or attractive (neg- 
ative dipole-dipole interaction energy) depending on the 
form of condensate cloud, its orientation relative to dipole 
polarization axes and trap geometry. Respectively sta- 
bility and collapse of BEC strongly depends on clouds 
anisotropy which opens a whole bunch of new phenomena 
to be observed and makes task of achieving and control 
of BEC especially challenging. 

Dipole-dipole interactions can dominate provided po- 
lar molecules are oriented by strong enough electric field. 
Similar effects can be achieved for ground-state atoms 
with electric dipole moments induced by a strong elec- 
tric field |(| 0], Another possible physical realization 
is atoms with laser induced electric dipole moments [||. 
Dipole-dipole interactions can be also essential in BEC of 
atomic gas with large magnetic dipole moments [Q, p"l| . 
Magnetic interactions are usually dominated by van der 
Waals forces but effects of magnetic interactions can es- 
sentially amplified by reducing of scattering length a via 
a Feshbach resonance jlOL 13 . Analysis of this Letter 



In this letter sufficient analytical criteria are developed 
both for catastrophic collapse of BEC of a trapped gas of 
dipolar particles and for long-time condensate existence. 
Sufficient criteria allows to predict condensate collapse 
or, opposite, its long-time existence for given condensate 
energy, E, number of particles, N, initial mean square 
width of condensate, and initial kinetic energy of con- 
densate. Analytical criteria are compared with results 
of variational approach ^] , where collapse was predicted 
based on the absence of local minimum of ground state 
of energy functional provided number of condensate par- 
ticle exceeds certain critical value. It is shown here that 
variational calculation gives threshold number of parti- 
cles and condensate energy which are located between 
parameters regions where analytical criteria predict col- 
lapse and long-time condensate existence, respectively. It 
is proven in this Letter that collapse certainly occurs pro- 
vided energy of the condensate exceeds a threshold value 
which is determined by the number of particles and trap 
parameters. Collapse of condensate is accompanied with 
dramatic contraction of the atomic cloud. Collapse is im- 
possible provided number of particles and initial kinetic 
energy of condensate are below the critical values. 

The time-dependent Gross-Pitaevskii equation (CPE) 
for atoms with long-range interactions and for a cylindri- 
cal harmonic trap is given by (6): 
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can be applied for both cases of electric and magnetic 
dipole-dipole interactions. 



where r = (x\ , X2, £3), VP is the condensate wave function, 
coupling constant g corresponds to short-range forces and 
is given by g — ATTh 2 a/m, a is the s-wave scattering 
length, m is the atomic mass, lu q is a trap frequency 
in X1X2 plane, and 7 is the anisotropy factor of the trap. 
*(r, t) is normalized to the total number of atoms in con- 
densate: N — J \^\ 2 d 3 r. It is assumed that the system 
is away from shape resonances of V(r) || and that the 
long-range potential is due to the dipole-dipole interac- 
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tion, and is given by 

[d 1 (r).d 2 (r')]-3[d 1 (r).u] [d 2 (r') • u] 



V(r-r') = 



(2) 



where u = (r — r')/|r — r'|. All the dipoles are assumed 
to point in the direction of trap axes (a^-direction), i.e., 
di = d 2 = da?3. Potential (||) with no dependence of d on 
r is a good approximation provided a typical interparticle 
distance exceeds a few Bohr radii. 

GPE (|l|) can be also obtained from variation of en- 
ergy functional, E : iti^ = j^z, where the condensate 
energy, 



E — Ek + Ep + Enl + Edd j 
is an integral of motion: = 0, and 



Ek= I ^\W\ 2 d 3 r, 
2m 



(3) 



E f 



mu>o(xf +X2 +7 x£)\i&fcrr, 



E 



NL 



\WdPr, 
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Consider time evolution of the mean square radius of 
the wave function, (r 2 ) = J r 2 \^\ 2 d 3 r/N. Using (|l|), in- 
tegrating by parts, and taking into account vanishing 
boundary conditions at infinity one gets for the first time 
derivative 



d t (r 2 ) = 
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J 2i Xj (^d Xj ^* - y*d Xj V)d 3 r, (5) 



where dt = Jj, d Xj = and repeated index j means 
summation over all space coordinates, j = 1, . . . ,3. 

In a similar way, after a second differentiation over t, 
one gets 
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Note that, in the case Ep = 0, V(v) = 0, Eq. (||) coin- 
cides with the so-called virial theorem for the GPE with 
local interactions Jl4|, [l5|, |l^, [l?], [l8| thus it is natural to 
call Eq. (H) by a virial theorem for GPE (|l|) . 

Using Eq. (j|) one gets (xjd Xj + Xjd x >.)V(r - r') = 
— 3U(r — r') and using Eq. (|J) one can rewrite virial 
theorem (||) as follows: 
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It is essential here that both local nonlinear term and 
nonlocal term are included into the energy E which is 
a conserved quantity. Catastrophic collapse of BEC oc- 
curs while (r 2 ) — * 0. From mathematical point of view 
it means that if, according to virial theorem (]?]), the 
positive-definite quantity (r 2 ) becomes negative in finite 
time then singularity in solution of Eq. ([!]) appears in 
a finite time before (r 2 ) becomes negative and singular- 
ity in solution of GPE occurs together with catastrophic 
squeezing of the distribution of \^\. Near singularity for- 
mation GPE is not applicable and another physical mech- 
anisms are important such as inelastic two- and three- 
body collisions which can cause a loss of atoms from the 
condensate ||. In addition, long term interactions are 
described by the dipole-dipole potential @ provided typ- 
ical distance between atoms in condensate exceeds a few 
Bohr radii. Note that regularization of potential (||) to 
avoid singularity at r = allows to prevent singularity 
formation in GPE |J, |(J. However GPE (Q) can still 
describe significant contraction of atomic cloud. 

Thus condition (r 2 ) — > provides a sufficient criterion 
of collapse of BEC. E.g. one immediately obtains from 
Eq. (0) that d 2 {r 2 ) < ^» and collapse is inevitable for 
E < 0. One can obtain however much more strict suffi- 
cient condition for collapse using generalized uncertainty 
relations between Ek , N, (r 2 ), dt(r 2 ) |ll| which follows 
from Cauchy-Schwarz inequality and Eq. (^|) with use of 
integration by parts = Re^, R= |*|): 
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Using Eqs. (Q), (||) one can obtain a basic differential 
inequality: 



d 2 (r 2 ) < 



1 



2mN 



12E 



Ti 2 / 97V m 2 N(d t {r 2 )) 2 



-10m^iVF(7)(r 2 } , (9) 



where F(-f) = 1 for 7 > 1 and F(j) = j 2 for 7 < 1. 
Change of variable, (r 2 ) = B A ' b /N gives the differential 
inequality: 



d 2 B < -- 
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which can be rewritten as 
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(11) 
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FIG. 1: Typical behaviour of potential U(B) from Eq. (|l|) 
for E < Ec-iticai (curve 1) E > E crit i Ca i (curve 2). Uo = 
(N 5 h 5 /m 5 uj ) 1/2 , B = (Nh/mujo) 5 ^. 



FIG. 2: Dependence of E c , var (curve 1) and E critica i (curve 
2) on trap aspect ratio 7 for N = N 0iVar . Both E C:Var and 
E cr iticai are given in units of Eq = h 7 ^ 2 Lo\j 2 /dPm?^ 2 . 



where 

U = -^-E&t* + ^?V 5 + f c Q 2 F (l)B\ (12) 
4m 3277T 8 

and f 2 (t) is some unknown nonnegative function of time. 
Equation (^lj) has a simple mechanical analogy jl6| with 
the motion of a "particle" with coordinate B under the 
influence of the potential force — dU Q^ in addition to the 
force —f 2 (t). Due to the influence of the nonpotential 
force —f 2 {t) the total energy £ of the "particle" is time 

B 2 

dependent: £(t) = -j- + [/(£?). Collapse certainly occurs 
if the "particle" reaches the origin B = 0. It is clear that 
if the particle were to reach the origin without the influ- 
ence of the force —f 2 (t) then it would reach the origin 
even faster under the additional influence of this nonpos- 
itive force. Therefore one can consider below the particle 
dynamics without the influence of the nonconservative 
force —f 2 {t) to prove sufficient collapse conditions. 

It follows from Eq. (|l2| ) that potential U (B) is a mono- 
tonic function for E < hu)oN[F( r y)5\ 1 / 2 /2 = E cri ti ca i 
(see curve 1 in Fig. 1) while for E > ^critical po~ 
tential U(B) has a barrier at B?,( 5 = 2,{e - [E 2 - 

1 /2 \ 

E l n ucai] )/[ 5muj o F h)] with particle energy £ m = 
U(B m ) at the top (see curve 2 in Fig. 1). One can 
separate sufficient collapse condition into three different 
cases: 

(a) for E < E cr iticai the particle reaches the origin in 
a finite time irrespective of the initial value of i?| Q ; 

(b) for E > Ecriticai and £(0) > £ m , the particle is 
able to overcome the barrier thus it always falls to the 
origin in a finite time irrespective of the initial value of 

(c) for E > Ecriticai and £ (0) < £ m , the particle is not 
able to overcome the barrier thus it falls to the origin in 
a finite time only if £?| t _ < B m . 

Note that it is proven here analytically only sufficient 
collapse conditions. It means that even if none of con- 
ditions a,b,c are satisfied one can not exclude collapse 



formation for some particular values of the initial condi- 
tions of Eq. ([!]). Generally it is determined by nonpo- 
tential force —f 2 (t). However inequality (^) reduces to 
equality for a Gaussian initial condition and 7 = 1: 

* = Arl/2^-3/4 (i 2 L3) -l/2 e -(^+^)/2^ e -^/2^ ^ 

in particular case with L3 = L p . 

One can compare sufficient collapse condition with 
results of Ref. B where collapse was predicted from 
variational analysis using Gaussian ansatz ( |l3"| ) to ap- 
proximate ground state of GPE for g — ([!]). It 
was concluded that collapse should occur provided en- 
ergy functional E has no local minima. A critical 
point was determined from the condition that local min- 
imum of energy functional E becomes a saddle point: 
(d 2 E/dL 2 ){8 2 E/dL 2 p ) = {d 2 E/dL 3 dL p ) 2 , 8E/dL 3 = 
d 2 E/dL p = 0. This allows to find critical number of 
particle, N c , V ar and critical value of energy functional, 
B c , var, as a function of system parameters L p , L3, 7, d.. 
Fig. 2 shows dependence of E c . V ar (curve 1) and E 'critical 
(curve 2) on trap aspect ratio 7 for N = N CtVar and 
g = 0. Note that the expression for E c<var , used in this 
Letter to draw curve 2 in Fig. 2, differs from Eq. (3) of 
Ref. Q. The authors of Ref. || already mentioned in 
the erratum pH ) that Eq. (3) of ||] is incorrect. How- 
ever a corrected formula was not given in the erratum 
pl[ . Explicit expression for E Ct var is not given in this 
Letter also because it is very bulky and will be given 
elsewhere. Gaussian ansatz ( |l3[ ) is not an exact solution 
of GPE (|l|) thus one can expect that actual critical value 
of energy E of ground state solution is lower. E cr iticai is 
determined here from sufficient collapse condition mean- 
ing that critical value of energy E of ground state solution 
is always above curve 2. One can conclude that actual 
critical value of energy is located between curves 1 and 
2. Accuracy of variation approximation can generally be 
obtained only from comparison with direct simulation of 
GPE (0). 

The Fourier transform of dipole-dipole interaction 
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FIG. 3: Energy lower bound E\ (in units of Eq) for TV > 
N c (curve 1) and N < N c (curve 2) versus Xf Xq . Xq — 
N 3 ^[mN Lo 2 /{^d 2 -3 5 )] 2/5 . 



potential (0) allows to find a sufficient condition of 
global existence (for arbitrary large time) of solu- 
tion of GPE (|l|). The dipole-dipole interaction en- 
ergy Ejjjj can be rewritten in k-space as Edd = 
(1/2) / | J R k | 2 l/ k d 3 k/(27r) 3 , where R k is a Fourier trans- 
form of |^| 2 and the Fourier transform of the dipole- 
dipole interaction in the limit of small atomic overlap 
distance is given by 0: Vk = — ^-d 2 (l — 3cos 2 a). 
Here a is the angle between k and d. Using inequal- 
ity Vic > -ifd 2 one gets E DD < -?fd 2 Y, where 
Y = J \V\ 4 d 3 r. Condition And 2 < ig results in E > 
for any particle number and collapse is impossible. Be- 
low it is assumed that And 2 > 3g. Y can be bounded as 



follows: Y < 337!^ N X I 2 X 3 I 2 , where X = J |VF| 2 d 3 r, 
and Nq = 18.94 is determined from a ground state so- 
lution, fa = XR(\r)e lX t , of nonlinear Schrodinger equa- 
tion: -\ 2 R + V 2 R + R 3 = 0, N = / R 2 d 3 r (see Ref. 
p3| ). Using these inequalities and Eqs. ©, (@), (§) one 
gets lower bound of energy functional: 



E > — X 
~ 2m 
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¥J 2 N~ 



N l/2 x 3/2 



For N > N c , N c = 2 3 / 2 3?i 5/2 Ay [5 5 / 4 (47rd 2 - 
3g)F(7) 1 / 4 m 3 / 2 w 1 / 2 ], the function E[(X) is a monotonic 
one (curve 1 in Fig. 3). For N < N c the function Ei(X) 
has a local minimum, E m i n (curve 2 in Fig. 3). Consider 
initial condition with 

N < N c , E mm <E< E max , Xi < X\t=o < X 2 , (15) 

where E max is a local maximum of Ei(X), and Xi,X 2 
are two of total three roots (X\ < X 2 < X3) of Eq. 
E = Ei(X). Any solution of GPE, corresponding to con- 
ditions ([l5]), will stay in the range X\ < X < X 2 at any 
time because regions below curves 1,2 in Fig. 3 are for- 
bidden for solution of GPE. One concludes that collapse 
is impossible in that case because collapse and singularity 
formation in GPE requires singularity in kinetic energy 



p2| , X — > 00. That could be understood e.g. from uncer- 
tainty relations (JsJ) . Eq. ( |l5| ) gives a sufficient condition 
of absence of collapse and in that case one can expect 
that energy functional E has a local minimum and sup- 
ports stable steady-state solutions. In original quantum 
mechanical problem that steady state is mctastable one 
because of finite probability of tunneling of condensate 
from local minimum which is outside the applicability of 
GPE and is not considered in this Letter. 

In conclusion, sufficient analytical criteria are devel- 
oped both for catastrophic collapse of BEC of gas with 
nonlocal long-range dipole-dipole interactions and for 
long-time collapse existence in the framework of GPE 
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